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ABSTRACT 

The last solar minimum has been unusually quiet compared to the previous minima (since space-based radiometric measurements 
are available). The Sun’s magnetic flux was substantially lower during this minimum. Some studies also show that the total solar 
irradiance during the minimum after cycle 23 may have dropped below the values known from the two minima prior to that. For 
chromospheric and coronal radiation, the situation is less clear-cut. The Sun’s 10.7 cm flux shows a decrease of ~ 4% during the solar 
minimum in 2008 compared to the previous minimum, but Ca n K does not. Here we consider additional wavelengths in the extreme 
ultraviolet (EUV), specifically transitions of Hei at 584.3 A and O v at 629.7 A, of which the CDS spectrometer aboard SOHO has 
been taking regular scans along the solar central meridian since 1996. We analysed this unique dataset to verify if and how the radiance 
distribution undergoes measurable variations between cycle minima. To achieve this aim we determined the radiance distribution of 
quiet areas around the Sun centre. Concentrating on the last two solar minima, we found out that there is very little variation in the 
radiance distribution of the chromospheric spectral line He i between these minima. The same analysis shows a modest, although 
significant, 4% variation in the radiance distribution of the transition region spectral line O v. These results are comparable to those 
obtained by earlier studies employing other spectral features, and they confirm that chromospheric indices display a small variation, 
whereas in the TR a more significant reduction of the brighter features is visible. 
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1. Introduction 

The minimum in solar activity between cycles 23 and 24 was dif¬ 
ferent in many aspects from other recent solar minima. It lasted 
longer, with nearly twice as many spotless days as during the 
previous eight minima, and was particularly deep and inactive 
(Schrijver et al. 2011; Tapping & Valdes 2011). The Sun’s in¬ 
terplanetary field in the recent minimum was at approximately 
half of the value it had during solar minimum 21 and was at its 
lowest level since the beginning of the measurements in 1963 
(Lockwood et al. 2009). This result is consistent with a halving 
of the Sun’s open flux (Vieira & Solanki 2010). 

A particularly interesting question is how the Sun’s radia¬ 
tive output reacted to the particularly low level of activity. Of 
the four main composites of total solar irradiance (TSI), ACRIM 
(Willson & Mordvinov 2003; Scafetta & Willson 2009), PMOD 
(Frohlich 2006; Frohlich 2009), RMIB (Mekaoui & Dewitte 
2008), and TIM (Kopp & Lean 2011), two (PMOD and the 
ACRIM) display a decrease in the TSI value during the mini¬ 
mum between cycles 23 and 24 compared to the previous min¬ 
imum (Woods 2010; Frohlich 2012). This decrease is indepen¬ 
dently obtained by models of solar irradiance (Yeo et al. 2014), 
which assume that the evolution of the surface magnetic field 
drives solar irradiance variability (Solanki et al. 2013). Thus for 
the TSI, which is formed mainly in the photosphere (Solanki & 
Unruh 1998), the evidence suggests a significant difference be¬ 
tween the two minima. 

The situation is less clear-cut for radiation formed in the up¬ 
per atmosphere. The solar radio microwave flux at 10.7 cm is 
commonly regarded as an activity indicator of the upper chro¬ 
mosphere and lower corona. The 10.7 cm radio flux has two 
main sources: thermal bremsstrahlung emission and gyroradi- 


ation; therefore, the enhancements in the temperature and the 
magnetic field of the Sun, as results from emerging sunspots, in¬ 
crease the 10.7 cm flux, thereby producing a high correlation be¬ 
tween the sunspot number and solar radio flux. However, in line 
with anomalous trends in some other solar activity indices during 
the last solar cycle, this well-pronounced correlation seems to 
have changed during solar cycle 23 (Tapping & Valdes 2011), es¬ 
pecially close to the solar minimum in 2008 (Svalgaard & Hud¬ 
son 2010). In addition, two different studies by Frohlich (2011) 
and Tripathy et al. (2010) show a very slight decrease (within 
4%) in the amplitude of the 10.7 cm flux over the last solar min¬ 
imum compared to the minimum in 1996. This is tiny compared 
with the changes displayed by other activity indicators between 
the two activity minima. 

Other solar activity indices, such as extreme ultraviolet 
(EUV) flux (Solomon et al. 2010; Woods 2010) and the Mg II 
index (Weber et al. 2013), also exhibit a similar decreasing trend 
over solar minimum 23. The measurements from the SOHO So¬ 
lar EUV Monitor (SEM) (Judge et al. 1998), for instance, reveal 
a 15% decrease with 6% uncertainty in the irradiance of the spec¬ 
tral band from 26 to 34 nm over minimum 23 with respect to the 
minimum 22 (Didkovsky et al. 2010). 

Particularly in the equatorial region, the solar surface is dom¬ 
inated by two main components: the quiet Sun and the active 
regions. The magnetic fields emerging from the solar photo¬ 
sphere through the boundaries between super granulation cells 
continually alter the upper solar atmosphere, resulting in a non- 
homogeneous EUV brightening (Fuller-Rowell et al. 2004). Al¬ 
though there are debates about various heating mechanisms of 
the chromosphere, a common view is that the main origin of 
the radiation from the chromosphere in the quiet Sun, espe- 
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cially in the internetwork regions, is the dissipation of the acous¬ 
tic waves, generated in the photosphere as a result of convec¬ 
tive motions of the plasma material at the top of the convec¬ 
tive zone (Rutten & Uitenbroek 1991; Carlsson & Stein 1997; 
Wedemeyer-Bohm et al. 2007; Schroder et al. 2012). However, 
the co-spatiality between chromospheric and photospheric net¬ 
work patterns, although not one-to-one (Schrijver et al. 1989; 
Loukitcheva et al. 2009), suggest that a magnetic heating mech¬ 
anism is at the root of the network formation in the higher atmo¬ 
sphere (Rutten 2012). 

The aim of this paper is to find out whether there is a change 
in the radiance of the upper chromosphere and transition re¬ 
gion of the quiet Sun network and internetwork between the ac¬ 
tivity minima following cycles 22 and 23. We present the re¬ 
sults of a statistical study of the radiance variability of the chro¬ 
mospheric He i 584.3 A (-20000 K) and transition region Ov 
629.7 A (-250000 K) spectral lines in the quiet Sun close to disc 
centre, over a period of 15 years from 1996 to 2011, with the 
main focus on comparing the last two solar minima. We used 
observations from the Coronal Diagnostic Spectrometer (CDS) 
aboard the Solar and Heliospheric Observatory (SOHO), which 
provides daily observations of the Sun in various wavelengths 
and has a long enough record covering the last two solar min¬ 
ima. To have a homogenous and consistent data set, all the mea¬ 
surements, during which an active region or bright feature was 
passing through the field of view of the instrument have been 
excluded from the final data set. This assures, to the maximum 
possible extent, that the results are based purely on observations 
of the quiet Sun. 

2. Observational data 

The CDS instrument (Harrison et al. 1995) aboard SOHO has 
been observing the Sun since 1996. The instrument consists of 
two spectrometers, the normal incidence spectrometer (NIS) and 
the grazing incidence spectrometer (GIS). The gratings in the 
NIS spectrometer disperse the incident radiation onto the NIS 
detector in two bands, covering the wavelength range from 310 A 
to 380 A (NIS-1) and from 517 A to 633 A (NIS-2). The spectral 
pixel size of the NIS ranges from 0.070 A at 310 A to 0.118 A at 
630 A. The spectral resolution was initially 0.35 A for NIS-1 and 
0.5 A for NIS-2 in terms of full width at half-maximum (FWHM) 
(Del Zanna et al. 2010). In June 1998 the contact with SOHO 
was lost due to operational errors. After recovery in Septem¬ 
ber 1998, the spectral profiles were broadened and asymmetric, 
resulting in a reduction in the spectral resolution of NIS (see 
Pauluhn et al. 2001). 

In this study, we used NIS daily synoptic observations (the 
so-called study program SYNOP_F). These observations cover 
the central meridian of the Sun from pole to pole. Each run of 
SYNOP_F consists of a series of nine raster scans, each cover¬ 
ing an area of 240" x 240" of the solar disc. Observations were 
taken employing the 2" x 240" slit (slit number 4) and by scan¬ 
ning along the solar v-direction in 120 separate exposures. The 
scale of the raster scans is 2" in solar v-direction and 1.68" in 
solar y-direction and each raster contains 120 x 143 spectra. In 
order to concentrate on the quiet Sun area close to disc centre, 
we only analysed the central rasters of each set of 9, with solar y 
pointings between -120" and +120". The complete data set we 
analysed covers a period of time from 1996 April 1 to 2011 May 
21 and consists of more than 4000 raster scans of daily measure¬ 
ments. As the main goal of this investigation is to look for any 
statistical variations in the radiance distribution of the quiet Sun 
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Fig. 1 . Appearance of a bright structure in the field of view of the CDS. 
The MDI magnetograms and EIT images were cropped in order to cover 
the same area of the solar surface as visible in the CDS images. Images 
acquired on 1996 July 27 


during the last two solar minima, we chose a time span of 1 year 
during each minimum and did a comparative study over various 
statistical parameters (as will be described in section 3). The 1- 
year periods were chosen to be from 1996 April 1 to 1997 April 
1 (minimum 22) and from 2008 July 1 to 2009 July 1 (minimum 
23). 

As mentioned above, a selection process was carried out to 
exclude the rasters scans containing any structure deemed to be 
unusually bright for the quiet Sun (e.g. due to decaying ARs). 
The selection was done in two ways: firstly, by visual inspec¬ 
tion of radiance images of CDS in the He i and O v lines and by 
referring to the He II 304 A images of the Extreme-Ultraviolet 
Imaging Telescope (EIT) (Delaboudiniere et al. 1995) and the 
Michelson Doppler Imager (MDI) (Scherrer et al. 1995) level 
1.8 full-disc magnetograms of the synoptic 96-minute series, and 
secondly, by analysing the goodness of the lognormal fits to the 
radiance distributions. This method will be explained in more 
details in section 3.3. Fig. 1 illustrates four images acquired by 
the three instruments on 1996 July 27. The two upper CDS im¬ 
ages show a small active region (AR) in the upper right corner. 
The same feature is also visible in the MDI and EIT images. Si¬ 
multaneous comparison of the four images makes it possible to 
identify the raster scans containing bright structures with high 
confidence. 

2.1. Data reduction 

The CDS raster scans were reduced using a set of SolarSoft 
routines. We applied VDS_CALIB for de-biasing, flat-fielding 
and conversion of instrument counts into physical units. We no¬ 
tice here that, in the case of spatially resolved observation, the 
measured physical quantity are spectral radiance (in units of 
mW sr _1 m -2 nm -1 ) and, after integration over wavelength, ra¬ 
diance (mW sr -1 m -2 ). However, it is customary (albeit formally 
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Fig. 2. Panel a is an intensity image in He i spectral line. Each vertical column in the image is corresponding to the average of the spectra over the 
wavelength dimension for a given x position (i.e. exposure). Panel b is a given exposure, corresponding to the white line in panel a, consisting of 
143 spectra along the y direction. Panel c is a spectrum recorded along the white line in Panel b. 


incorrect) to use the word intensity in the place of radiance. We 
will also make use of the term intensity, always meaning radi¬ 
ance. The cosmic ray hits were removed by applying the stan¬ 
dard routine CDS_NEW_SPIKE. The CDS\NIS data are stored 
in a cubic structure with the dimensions of x, y and A, corre¬ 
sponding to solar x and y directions and the wavelength, meaning 
that each spatial "pixel" of a given raster scan contains a spec¬ 
trum (see Fig. 2). We analysed individual spectra and excluded 
the ones which are highly damaged due to cosmic ray hits. On 
average, in a typical raster scan, approximately 3% of the spec¬ 
tra were flagged out. NIS spectra also suffer from geometrical 
distortions such as rotation and tilt of the spectra, resulting from 
the misalignments between the detector and the grating, and be¬ 
tween the grating and the slit, respectively. These effects were 
corrected by applying the NIS_ROTATE routine. After apply¬ 
ing the instrumental standard calibration and correction routines, 
the solar radiance at each spatial pixel was determined by least- 
squares-fitting of a single Gaussian function with constant back¬ 
ground to a pre-loss spectrum, i.e. any spectrum recorded prior 
to June 24, 1998, and by fitting a broadened Gaussian line profile 
to each post-loss spectrum 1 and finally by extracting the fitting 
parameters. The final output of the whole process is a set of daily 
intensity images, with the intensity at each pixel corresponding 
to the area under the Gaussian after subtracting the background. 

2.2. CDS point-spread function and data binning 

By intercalibration of SUMER detector A and CDS\NIS, 
Pauluhn et al. (1999) pointed out that CDS has a broad PSF cor¬ 
responding to FWHM X = 6" and FWHM y = 8", which is almost 
three times larger than the pixel size of the NIS detector in the 
x direction and over four times larger in the y direction. This 
study was based on observations prior to SOHO loss of control 
and the post-recovery data are believed to have an even larger 
point-spread function. Given the broad PSF, binning the inten¬ 
sity images sounds logical, since it increases the signal to noise 
ratio without reducing the spatial resolution. We binned NIS data 

1 See CDS Software Note Nr.53: http: //solar.bnsc. rl. ac. uk/ 
swnotes/cds_swnote_53.pdf 
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Fig. 3. Illustration of long-term drop in CDS\NIS responsivity. The 
CCD readout bias has been removed from the data right after reading 
them in. No other calibration has been applied to the data. The solid 
lines and their associated error bars represent the mean and the standard 
deviation of the average intensity values over 90-day time intervals. The 
vertical dashed-lines indicate the two one-year time spans which we 
used to compare the two solar activity minima. 


with six different binsize values: from binsize equal to 2 (binning 
over groups of 2 x 2 pixels) up to binsize equal to 7 (binning over 
groups of 7 x 7 pixels). The down-sampling (binning) of the data 
was carried out after applying instrument calibration corrections 
and prior to spectra fitting. In this way, we prepared seven data 
sets corresponding to the different binsize values. The binsize 1 
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Fig. 4. Long-term variation in daily average intensities of the He i and 
O v spectral lines. The solid lines and their associated error bars repre¬ 
sent the mean and the standard deviation of the average intensity values 
over 90-day time intervals before (grey) and after (black) applying the 
Del Zanna et al. correction factors. The vertical dashed-lines indicate 
the two one-year time spans chosen to cover the solar activity minima. 

in the plots in this paper, means that the images have not been 
down-sampled and are in the original format. 

3. Analyses and results 

3.1. Daily average intensities 

The CDS/NIS suffers from long-term degradation and decrease 
in responsivity. To verify the drop in the instrument responsivity, 
we read in the raw data, removed CCD read-out bias and be¬ 
fore applying any further calibration, averaged them over spatial 
and spectral dimensions. Panel a in Fig. 2 displays an intensity 
image in the He i spectral line in debiased_ADC unit, i.e. after 
debiasing and prior to applying calibrations. Each x position on 
the image corresponds to the average of an exposure over the A 
dimension. Panel b shows a given exposure corresponding to v 
position number 30 of the intensity image (white vertical line in 
Panel a). Panel c depicts a spectrum recorded along the white 
line in panel b at y position number 101, meaning that the value 
given to the pixel coordinate [30,101] in panel a is the average 
of the spectrum in panel c. The daily average intensities with the 
unit of "debiased-ADC" were calculated by averaging over x, y, 
and A dimensions. These values are illustrated in Fig. 3, reveal¬ 
ing a drastic drop in the responsivity of the instrument over the 
years. To compensate for this systematic drop, a set of calibra¬ 
tion and correction routines has been produced and updated by 
the CDS team. 

After applying standard calibrations to the data, the daily av¬ 
erage intensities in physical units were calculated by fitting the 
spectra (as explained in section 2.1) and averaging the intensity 
in the spectral line over the entire intensity images. The grey line 
in Fig. 4 illustrates the long-term variations in average intensities 
in the He i and O v spectral lines after applying standard calibra¬ 


tions. It represents the mean intensity values over 90-day time 
intervals and the vertical bars represent its corresponding stan¬ 
dard deviation, respectively. Del Zanna et al. (2010) suggested 
that the instrument standard calibration routines overcorrect the 
long-term drop in the responsivity of the instrument, resulting 
in a steady increase in the calculated radiance. By analysing the 
NIS daily synoptic observations from 1996 to 2009, they intro¬ 
duced a new set of correction factors to fix this issue. The black 
solid line and its corresponding error bars in Fig. 4 display the 
90-day averaged intensities after applying the Del Zanna et al. 
(2010) correction factors. 

The average Hei intensity values calculated by using only 
the standard corrections (grey line in Fig. 4), demonstrate an in¬ 
crease starting approximately from 2005. Such a trend is not seen 
in O v. Applying Del Zanna et al. correction factors (black lines) 
moderates the issue to a great extent for He i but then the O v in¬ 
tensity values seem to be overcorrected for the period after 2008. 

The correction factors introduced by Del Zanna et al. have 
been determined based on the assumption that the variations in 
the radiances of the chromospheric and transition region spec¬ 
tral lines are mainly due to the existence of the ARs and that 
the mean radiances of the spectral lines in the quiet Sun at disc 
centre essentially does not respond to the solar activity cycle. Al¬ 
though it is expected the solar radiance to be of the same order of 
magnitude at different solar minima, this assumption essentially 
precludes any investigation such as the present one, which aims 
at detecting even small variations in the quiet Sun’s radiance, as 
long as they are significant. The small differences in radiance 
properties of the quiet Sun between two solar minima can be re¬ 
vealed only if the solar activity parameter under investigation is 
adequately independent of instrumental degradation. The aver¬ 
age intensity values do not seem to satisfy this pre-requisite. The 
alternative is to identify and study parameters of the intensity 
distribution which are not affected by instrumental degradation. 
This is the approach we decided to follow. 

3.2. Contrast ratios 

In order to minimize the effects of the calibration factors in the 
analysis, we introduce two ratios, the 5 and 10 percent contrast 
ratios (called 5% and 10% CR in the following), defined as the 
ratios between the average intensities of the brightest 5% and 
10% pixels within an intensity map and the average intensities 
of the darkest 5% and 10% pixels, respectively. Fig 5 displays a 
histogram corresponding to a given intensity map. The brightest 
and darkest pixels are highlighted in red and blue, respectively. 
The advantage of defining these CRs is that the correction factors 
which have already been applied to the data will be cancelled out 
by dividing. From now on we confine our analysis to the solar 
minima. Fig. 6 illustrates the CRs for three different binsizes (see 
section 2.2). In each panel of the figure the four colours distin¬ 
guish between CRs of the two spectral lines at each minimum. 
The x axes of the bottom and top of each column of the figure 
represent the date during solar minima 22 and 23, respectively. 
Excluded from the final plotted dataset, are days on which the 
CDS images sampled not purely quiet Sun (see Fig. 1) and/or 
because substantial parts of the image were missing. 

Obviously the 5% CRs are always larger than the 10% CRs. 
This is not surprising, since averaging over a larger number of 
pixels (10% vs. 5%) at the bright and dark extremes of the inten¬ 
sity distribution results in closer mean values and consequently a 
closer-to-unity ratio. In addition, the CRs decrease as the binsize 
increases. This is also expected due to the fact that the inten¬ 
sity images tend to reduce towards grey images as more pixels 
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Fig. 5. Histogram gives the number of pixels displaying intensity within 
a certain range on 1996 June 6. The 5% pixels with the lowest intensity 
are shaded dark blue. These plus the lightly shaded blue pixels represent 
the 10% darkest pixels. Similarly the more strongly shaded red pixels 
represent the 5% brightest pixels. The 10% brightest pixels are the sum 
of the pixels shaded in dark and light red. 


are binned together. It is also evident that for all binsizes O v 
shows higher CRs than He i. This difference starts from approxi¬ 
mately 50% higher CR values for binsize 1, and drops as binsize 
increases. The possible reasons for these differences will be dis¬ 
cussed in detail in the discussion section. What is more in focus 
here is the differences in CRs of the same spectral line over the 
two solar minima, e.g. differences between the solid and dotted 
horizontal lines (for O v) or between the dashed and dash-dotted 
horizontal lines (for Hei). For a binsize of 1, this difference is 
quite large for both spectral lines. However with increasing bin- 
sizes, it diminishes in He i, but persists in O v. 

To have a direct comparison between the two minima, we 
computed the ratio of the mean CR of minimum 23 to that of 
minimum 22 ( mm23 CR / min22 CR) for both spectral lines. This ra¬ 
tio of CRs is illustrated in Fig. 7 as a function of binsize. The 
ratios are either close to or below the unity-line for all binsizes 
except binsize 1. This is related to the fact that the darker pixels 
are dominated by noise for the data taken during solar minimum 
23 due to the drop in the instrument responsivity (see Fig. 3). 
Therefore, when making Gaussian fits to individual spectra in 
darker areas, the signal to noise ratio is so low that often what 
is being fitted is noise rather than the spectrum. Therefore, the 
line intensities in these areas tend to be underestimated. This ef¬ 
fect is larger for the data of minimum 23 and explains why the 
CRs are larger during solar minimum 23. In addition to binsize 
values, the plots in the figure also highlight the effect of data ex¬ 
clusion. Three different versions of data exclusion were applied 
to examine how sensitive the ratios are to the data exclusion. In 
version 1, almost 10% of the images containing the most obvious 
bright structures (chosen visually) were excluded from the data 
set. Versions 2 and 3 correspond to 20% and 30% data exclu¬ 
sion, respectively. They represent progressively stricter rejection 
criteria. It is evident from the plots that by excluding 20% of 


the data, the He i has already reached a stable state since the ra¬ 
tios for versions 2 and 3 lines almost overlap and excluding more 
data points does not change the minimum-to-minimum ratios for 
any binsize. This is different for O v, where excluding more data 
keeps pushing the ratios toward the unity line. 

The plots also show that the minimum-to-minimum ratios are 
less than 1 for all binsizes (except binsize 1) in case of O v and 
are very close to unity in He i. This suggests that the intensity 
distribution of O v underwent significant changes between the 
last two solar minima while that of the upper chromospheric He i 
did not. 

3.3. Intensity distribution of the quiet Sun radiation and 
lognormal fitting 

Various studies have shown that the intensity distribution of 
the quiet Sun in the chromosphere and transition region has a 
skewed shape, with a pronounced peak at relatively low intensi¬ 
ties and an extended tail at higher intensities. This skewed distri¬ 
bution has been modelled mainly by either a sum of two Gaus¬ 
sian functions (Reeves 1976, Gallagher et al. 1998) or a single 
Lognormal function (Griffiths et al. 1999, Pauluhn et al. 2000, 
Fontenla et al. 2007), each representing a different hypothesis 
about the source of the radiation in quiet Sun regions. These 
two models imply different heating scenarios. A two-component 
function assumes two intrinsically different heating mechanisms 
for radiation originating in the solar network and internetwork, 
while a single-component model assumes similar heating mech¬ 
anisms for both network and internetwork (Pauluhn et al. 2000). 
The common agreement today is that the quiet Sun radiation fol¬ 
lows in general a lognormal distribution function. Thus Pauluhn 
et al. (2000) showed that despite its fewer free parameters, a log¬ 
normal reproduced the observed distributions significantly better 
than the sum of two Gaussians. 

It is believed that a system variable can be modelled with a 
lognormal distribution if it is a multiplicative product of many 
independent forces acting on that system (Gronholm & An- 
nila 2007, Limpert et al. 2001). These expressions can be ap¬ 
plied to solar physics through, e.g., the nanoflare heating model, 
which is a potential mechanism for coronal heating (Parker 1983, 
1988). The consistency between lognormal distribution of tran¬ 
sition region radiation at the quiet Sun and nanoflare heating of 
the corona has been studied by Pauluhn & Solanki (2007) and 
Bazarghan et al. (2008). 

As mentioned earlier, in order to study the variation of the 
radiance distribution over time, we fit the individual distribu¬ 
tions with a lognormal function, extract the best fit parameters 
and investigate their long-term variability. The lognormal distri¬ 
bution function can be mathematically defined by replacing the 
independent variable (here intensity I) in the normal distribution 
function with its logarithm. In other words, if log (/) is normally 
distributed, the variable I is lognormally distributed: 


P(I) = 



(In (/)-yu) 2 \ 
2cr 2 )' 


( 1 ) 


The above equation shows a lognormal probability den¬ 
sity P(I) as a function of intensity. Here yu = (In (/)), cr = 
VVar(ln (/)) and Vo is the normalization factor. The fi and cr are 
the mean and standard deviation of ln(7). We use natural loga¬ 
rithm in the definition of the lognormal function. Using another 
base to the logarithm only rescales the fitting parameters, but 
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Fig. 6. 5% (left column) and 10% (right col¬ 
umn) contrast ratios (CR) vs. time. See the main 
text for the definition and Fig 5 for an illus¬ 
tration. The one-year periods during minimum 
22 and minimum 23 are given by the bottom 
and top axis, respectively. Each panel repre¬ 
sents results for a different binsize (labelled in 
the panel). The horizontal lines show the aver¬ 
age value of each CR. Binsize 1 means that the 
intensity images are not binned. 
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Fig. 7. Ratio between CR at the two solar minima for different binsizes. 
The blue horizontal dotted line indicates the unity line. The dashed, 
dotted and solid lines are corresponding to the three versions of data 
exclusion, corresponding to 10%, 20% and 30% data exclusion due to 
existence of bright features inside the field of view, according to pro¬ 
gressively stricter criteria. The error bars represent the lcr error in the 
ratios. 

eventually returns the same family of distributions. The back- 
transformed parameters are then defined as 

jf = , cr* = e (T (2) 

and called scale parameter and multiplicative standard devi¬ 
ation or shape parameter, respectively (Limpert et al. 2001). 

Fig. 8 illustrates two intensity images in the O v line, their 
corresponding histograms and two MDI magnetograms of a 
slightly larger solar area, encompassing the area scanned in O v. 
The upper row corresponds to a date on which part of an ac¬ 
tive region passed through the field of view, while the bottom 
row depicts quiet Sun conditions. Each plot consists of two his¬ 
tograms: the actual intensity histogram shown by the solid black 


line and the normalized histogram of the natural logarithm of the 
intensity (log(/)) represented by the solid grey line. The dashed 
lines are Gaussian fits to the normalized histogram of the log¬ 
arithm of the intensity distribution. A comparison between the 
two plots reveals how significantly the distributions can deviate 
from a lognormal distribution when a particularly bright struc¬ 
ture, such as a part of an AR, appears in the field of view, x 2 
values, as will be explained later, were applied to exclude the 
non-quiet Sun data from the analyses as an alternative approach 
to visual inspection. 

The scale and the shape parameters of the individual inten¬ 
sity distributions were calculated according to the p and cr values 
of the normal fits, as formulated before. Note that the scale pa¬ 
rameter is sensitive to the instrumental corrections. Referring to 
the mathematical definitions of the shape and scale parameters 
given above, multiplying an intensity image to a given number 
(as a correction factor) shifts the position of the peak of the his¬ 
togram and consequently changes the scale parameter, but not 
the shape parameter. This can be seen in Fig. 9, which exhibits 
the variation in scale and shape parameters of the intensity dis¬ 
tributions of the He i spectral line. The scale parameter shows a 
similar trend to that of the average intensities of He i shown in 
Fig. 4, while the shape parameter does not depict any systematic 
trend. 

The plot of the shape parameter in Fig. 9 also demonstrates 
an increase as we approach the solar maximum. This happens 
since the overall emergence of active areas during the solar max¬ 
imum gives rise to the number of bright pixels in the tails of 
the radiance distributions. Schiihle et al. 2000 found similar in¬ 
creasing trends in the radiance of various ultraviolet lines be¬ 
tween 1996 and 2000. They showed that the network and inter¬ 
network radiances increase uniformly over the solar cycle. How¬ 
ever, whether such an increase occurs indeed uniformly over net¬ 
work and internetwork can only be seen in daily average intensi¬ 
ties, such as those in Fig. 4 (from the minimum to the maximum 
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Fig. 8. From left to right: MDI magnetograms, O v integrated intensity scans and intensity distribution functions obtained from the intensity 
images. The solid black lines in panels c and f show the intensity histograms on a linear intensity scale scale (corresponding to the axes on bottom 
and left side of panels c and f). The solid red lines are the normalized histograms of the logarithm of the intensity (corresponding to the axes on 
top and right side of panels c and f). The dashed lines are lognormal fits to the normalized histograms. The surface area visible in CDS images is 
marked by the white rectangles in the MDI magnetograms. 


of cycle 22). However, we can not verify this finding due to the 
uncertainties in the absolute intensity values. 

Various statistical parameters can be defined using the /i and 
the cr values. We do not analyse the parameters that are depen¬ 
dent on id (e.g., scale parameter, mean and mode) for the same 
reason as we don’t analyse the average intensities. These param¬ 
eters can be used to improve the instrument calibration, as was 
done by Del Zanna et al. (2010). We decided to follow a differ¬ 
ent approach here. Parameters that only depend on cr are best 
for studying long-term variations. Four such statistical parame¬ 
ters were chosen for further study, the shape parameter cr *, the 
skewness y, the coefficient of variation (CV) and the ^ 2 value of 
the fits. The CV parameter is a measure of the variable disper¬ 
sion, while the skewness is a measure of the asymmetry of the 
distribution. In general, the distribution is more asymmetric and 
dispersed for larger cr, while a smaller cr implies that the distri¬ 
bution is more symmetric. Skewness and CV are computed as 
follows: 

y ~ ( e 0-2 + 2) V ^ 2 - 1 (Skewness) (3) 

CV = V ^ 2 - 1 (Coefficient of Variation) (4) 

Fig. 10 displays the ratios of the minimum of cycle 23 values 
to those during cycle 22 minimum of various parameters for dif¬ 


ferent binsizes. Plotted are such ratios of cr*, 7 and CV averaged 
over the periods of time defined as minima. Excluding the x 2 
plot which gives a measure of the goodness of the fits, all other 
plots depict lower than unity ratios for all binsizes greater than 1 , 
specially for O v. This suggests that the intensities follows a less 
asymmetric distribution during solar minimum 23, i.e., the inten¬ 
sity distributions are closer to a normal distribution. This effect 
is less pronounced in He 1 , with closer-to-unity ratios compared 
to the O v spectral line. This implies that the mean values at the 
two ends of the intensity distributions are closer in magnitude, so 
that this result is in agreement with the 5% and 10% CRs being 
smaller at minimum 23, as shown in Fig. 7. Here again the ex¬ 
ception is the binsize 1 , which shows a similar behaviour to that 
of the CRs. The^f 2 plot helps to determine which amount of bin¬ 
ning is most appropriate for comparing the two minima (where 
thex 2 is the same in both minima). The best consistency between 
the data for two minima, in terms of signal to noise ratio, can be 
achieved for binsizes 3 and 4 for both spectral lines. 

As already shown in Fig. 8 , the intensity distributions of the 
images with any trace of activity exceeding that of the typical 
quiet Sun deviate from a lognormal. The ^ 2 values of the lognor¬ 
mal fits, as mentioned before, can therefore also be utilized as a 
potential reference for excluding scans that do not cover purely 
quiet Sun. Fig. 11 shows the minimum-to-minimum ratios of the 
5% and 10% CRs, similar to the one which is already shown 
in Fig. 7, but this time instead of excluding the data by visual 
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Fig. 9. Long-term variations in the scale (upper panel) and shape (lower 
panel) parameters of the lognormal fit to the intensity distribution of the 
He i spectral line. The solid lines and vertical error bars correspond to 
90-day average values and standard deviations, respectively. 


inspection, we selected them according to the^ 2 of their lognor¬ 
mal fits. Here again three different versions were introduced cor¬ 
responding to the number of data points which were excluded. 
0,10, 20 and 30% in the images label implies that this fraction of 
the intensity images with highest^ 2 values, was removed before 
the CRs were computed. In order to check if the ratios are sen¬ 
sitive to the histogram bin width and/or to the weighting of the 
bins for lognormal fitting, we tested two different methods for 
defining the histogram bin width and two methods to weight the 
histogram bins. The bin widths were either given a fixed value 
equal to 20 mW/m 2 /sr for the entire data set and for both spec¬ 
tral lines or alternatively determined according to Scott’s normal 
reference rule for each individual histogram. Scott’s normal ref¬ 
erence rule estimates the optimal bin width for a histogram by 
assuming the underlying distribution to be a Gaussian and by 
minimising the integrated mean square error (IMSE) of the his¬ 
togram (Scott 1979). For lognormal fitting, the bins were either 
given a constant weight or a statistical (Poisson) weight. The 
four different combinations of these methods for the two CRs 
are illustrated in Fig. 11. In most cases, there is a remarkable 
gap between the 0% line and the other lines, whereby ’0% line’ 
refers to the ratio of CRs obtained without removing any data 
(i.e. keeping also the brightest scans, most strongly affected by 
active regions). The fact that the 10, 20 and 30% lines almost 
overlap each other, implies that the ratio between the CRs ap¬ 
proaches a saturation level when going to larger exclusion per¬ 
centages. This was less evident for Ov in Fig. 7, where the ex¬ 
clusion was done by visual inspection. Here again, regardless of 
which method for determining the bin widths and bin weights is 
used, the He i line does not show any significant difference be¬ 
tween the two minima. The O v CRs, on the other hand, exhibit a 
drop of approximately 8% during minimum 23 compare to that 
of minimum 22, which is almost twice the change we saw in the 
CRs when excluding the bright days by visual inspection. 


4. Discussion and conclusion 

The degradation of the CDS instrument’s photometric sensitivity 
makes the absolute intensity of spectral lines unreliable for long¬ 
term investigations of subtle solar changes such as the difference 
in brightness between the minima following cycles 22 and 23. A 
qualitative study of the parameters independent of absolute in¬ 
tensity, is a possible alternative that we consider here. For this 
purpose, we defined two contrast ratios (CRs), as the ratio be¬ 
tween the mean intensity of the brightest and darkest parts of 
the intensity images (see Fig. 5) . This way we cancelled out the 
effect of correction and calibration factors. We also fitted log¬ 
normal distribution functions to the intensity distributions and 
extracted the fit parameters. Eventually a comparative study of 
all these parameters over the last two solar minima was carried 
out to investigate the variations in the radiance distribution. We 
also excluded the images with bright structures in the field of 
view to concentrate only on the quiet Sun area. The exclusion 
procedure was based either on visual inspection (Fig. 1) or on 
the goodness of the lognormal fits to the intensity histograms 
(Fig. 8). In addition, the NIS raster scans were down-sampled 
(binned) prior to any further analyses to compensate for the ef¬ 
fect of reduced sensitivity of the instrument by increasing the 
signal to noise ratio over darker areas. This does not effectively 
reduce the instrument spatial resolution up to binsize 3, as the 
instrument PSF is almost three times larger than the CDS actual 
pixel size (see section 2.2). Therefore, binsizes 3 and 4 are be¬ 
lieved to be the optimal values for down-sampling in order to 
have maximum signal to noise ratio and minimum loss of spatial 
resolution. This is in agreement with previous work by Brkovic 
et al. (2000). 

The CRs in O v exhibit a remarkable drop over solar mini¬ 
mum 23. The magnitude of the reduction, however, is sensitive to 
the data exclusion methods. From Fig. 7, when considering bin- 
sizes 3 and 4, exclusion of 30% of the images by visual inspec¬ 
tion reveals a drop of around 4% in the minimum-to-minimum 
trend. The drop is almost doubled when the exclusion is based 
on the^ 2 of the lognormal fits, displayed in Fig. 11. The latter 
figure shows that no substantial changes occur when rejection 
levels larger than 10% are used. 

The CRs, are in principle a measure of contrast between the 
network and internetwork radiances. The increase of the ratio of 
the CRs as the number of excluded bright images is increased, is 
due to the fact that we are left with progressively "purer" quiet 
Sun when more and more images with bright features are ex¬ 
cluded. The Ov radiance distributions undergo a shift from a 
skewed form to a more symmetric shape during the very quiet 
minimum of cycle 23. Such a variation of the distribution man¬ 
ifests itself in terms of a decrease in the shape and skewness 
parameters illustrated in Fig. 10, where the skewness of the in¬ 
tensity distribution of the O v line shows a drop of 3.5% over so¬ 
lar activity minimum 23. The smaller dispersion in O v intensity 
distribution, as quantified by the coefficient of variation (panel c 
of Fig. 10), during activity minimum 23 suggests that either the 
thermal structure of the internetwork has changed between the 
minima following cycles 22 and 23, or the area covered by the 
network or the brightness of the network. The latter two possible 
explanations relate to the amount of magnetic flux. Lower mag¬ 
netic flux is also found to be the cause of the lower TSI (we refer 
to TSI composites PMOD and ACRIM, although the RMIB and 
TIM composites do not show a decrease) during the minimum 
of cycle 23 (cf. Ball et al. 2012; Yeo et al. 2014). On the other 
hand, the studies carried out by Frohlich (2011) and Tripathy 
et al. (2010), show also a drop in the 10.7 cm flux up to 4% over 
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Fig. 10. Cycle 23 activity minimum to cycle 
22 activity minimum ratios of various statistical 
parameters of the lognormal fits. The solid and 
dash-dotted lines represent the ratios in O v and 
He i spectral lines, respectively. The error bars 
represent the lcr error in the ratios. 


minimum 23 with respect to the previous one. We note that the 
10.7 cm radio flux is a global solar activity factor and it is not 
trivial to compare it directly with our results, which are entirely 
based on local observations of the solar disc centre. 

The higher CRs in transition region Ov line during solar 
minimum 22 can also be attributed to the presence of active net¬ 
work elements, being dispersed within the quiet Sun. The ac¬ 
tive network results from the decay of active regions (Skumanich 
et al. 1984; Worden et al. 1998, 1999) and may last for several 
solar rotations before diffusing into the quiet Sun (Woods et al. 
2000). The exclusion criteria for rejecting images with bright 
structures may have caused us to miss the active network contri¬ 
bution in the quiet Sun radiance distribution. 

The He i spectral line, in contrast to the O v line, shows al¬ 
most no minimum-to-minimum variations. This is evident when 
referring to the CRs in Figs. 11 and 7. Both 5% and 10% CRs 
at binsizes 3 and 4 are very similar in the last two solar min¬ 
ima and their minimum-to-minimum ratios are within 1 % of the 
unity line. It is also interesting to note that the He i line is less 
sensitive to the exclusion of data due to the presence of magnetic 
activity. Excluding only 10% of the images with the most obvi¬ 
ous bright structures, is sufficient for the total brightness not to 
be affected by the remaining traces of activity. This is because 
the contrast between network and internetwork is larger in the 
transition region O v line than in the chromespheric He i line. In 
addition, Woods et al. (2000) show that the contrast for the active 
network is more enhanced for transition region emission (such as 
in He n 304 A) than it is for chromospheric or coronal emission. 
This is also the likely reason why the O v intensity distribution 
has a more skewed shape over minimum 22. 


These findings are compatible with those of Livingston et al. 
(2010), who found that the chromospheric Ca II K line, mea¬ 
sured over solar disc centre, is constant (within 1%) between 
the last three solar activity minima. This is important, for the 
Ca II index is usually regarded to display a relatively tight rela¬ 
tionship with the photospheric magnetic field (Skumanich et al. 
1975; Schrijver et al. 1989; Lukianova 2009). The He i and Ca II 
K results imply a certain level of decoupling between the chro¬ 
mospheric radiation and the photospheric field at least at very 
low level of activity. 

In a recent study by Andretta & Del Zanna (2014), the au¬ 
thors also analysed the radiance distribution over solar cycle 23 
by investigating the full-disc and central meridian observations 
of the CDS studies USUN and SYNOP, respectively. They con¬ 
cluded that the quiet Sun radiance has undergone no significant 
changes during solar cycle 23, regardless of which spectral line 
is considered. Although we came to the same conclusion for He i, 
we found a minimum-to-minimum variation in the radiance dis¬ 
tribution of the O v. Nevertheless, the values they found as the 
average widths of the lognormal fits over the full cycle are in 
good agreement with our findings (when rescaled from natural 
logarithm to normal logarithm). We measured average width val¬ 
ues for histograms in Hei and O v to be equal to 0.18 ± 0.01 and 
0.22 ± 0.02 dex, respectively. The values reported in the study of 
Andretta & Del Zanna (2014) are 0.16 ± 0.015 dex for Hei and 
0.19 ± 0.02 dex for Ov. The small difference might be due to 
the specific fitting approach applied by (Andretta & Del Zanna 
2014), in which only the peak of the histograms, corresponding 
to the internetwork area, are fitted and the contribution of the 
brighter network is ignored. 
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Fig. 11. Minimum-to-minimum ratio of the 5% 
and 10% CRs. Each panel, consists of the ratios 
for both spectral lines and for different amounts 
of data exclusion, based on the x 2 values of a 
lognormal fit to radiance distributions. The ef¬ 
fect of two different methods of determining 
the histogram binsize and weighting them for 
fitting has been emphasized. The abbreviations 
cut and set and pois are for constant, Scott’s 
rule and Poisson, respectively. The differently 
coloured lines are labelled in panel b. 


In the model introduced by Dowdy et al. 1986, a two- 
component mechanism is assumed to be responsible for pro¬ 
ducing the emission in the transition region (TR): a downflow 
of the heat within the magnetic funnels which expand upwards 
throughout the network structures connecting to corona, and the 
internal heating of the gas inside low-lying small loops, which 
are magnetically insulated from the corona. Each of these two 
components connects the radiance formation in the TR to the 
chromosphere and corona. 

The difference in behaviour of the chromospheric and TR 
lines found here may have different origins. One possibility is 
that the He i line reacts less strongly to the presence of a mag¬ 
netic field than the O v line. Another is that the heating of the 
quiet chromosphere is much more strongly determined by acous¬ 


tic waves than the TR. As a result, in He i we are basically seeing 
basal flux in the very quiet Sun, while in TR lines a more magnet¬ 
ically dominated heating mechanism is active. The former does 
not change from minimum to minimum, while the latter does. 
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